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Abstract. Sulfur abundances have been determined for ten stars to resolve a debate in the literature on the 
Galactic chemical evolution of sulfur in the halo phase of the Milky Way. Our analysis is based on observations 
of the Si lines at 9212.9, 9228.1, and 9237.5 A for stars for which the S abundance was obtained previously from 
much weaker Si lines at 8694.0 and 8694.6 A. In contrast to the previous results showing [S/Fe] to rise steadily 
with decreasing [Fe/H], our results show that [S/Fe] is approximately constant for metal-poor stars ([Fe/H];$— 1) 
at [S/Fe] ~ -1-0.3. Thus, sulfur behaves in a similar way to the other a elements, with an approximately constant 
[S/Fe] for metallicities lower than [Fe/H] ~ — 1. We suggest that the reason for the earlier claims of a rise of [S/Fe] 
is partly due to the use of the weak S I 8694.0 and 8694.6 A lines and partly uncertainties in the determination 
of the metallicity when using Fel lines. The Si 9212.9, 9228.1, and 9237.5 A lines are preferred for an abundance 
analysis of sulfur for metal-poor stars. 



1. Introduction 

Galactic chemical evolution of the a elements (O, Ne, Mg, 
Si, S, Ar, and Ca) might be expected to be similar from 
element to element. It is thought that their production in 
the early Galaxy was dominated by synthesis in massive 
stars and ejection into the interstellar medium by super- 
novae (SN Type II). Then, relative abundances of a ele- 
ments and of their abundances with respect to iron reveal 
information about the yields from SN II and the depen- 
dence of those yields on the initial metallicity. After a 
passage of time, the Galactic gas was contaminated by 
ejecta of Type la supernovae (exploding white dwarfs). 
These produce little of the a elements but large amounts 
of iron-group nuclei. Thus, the relative abundance of a 
elements to iron declined. 

Accurate data on the abundances of a elements for 
stars in the halo and disk are essential if one is to 
extract information about the nucleosynthesis of these 
elements. First estimat es of the S abunda n ce in halo 
stars were provided by pTakada-Hidai et al ] ll2002l) and 
llsraelian fc Rebolol l)200lj) who found that [S /Fe] increased 
approximately linearly with decreasing [Fe/H]: [S/Fe] ~ 
-0.4 X [Fe/H] from the Sun ([Fe/H] = 0) to [Fe/H] 
-3 , the m etallicity limit of the samples. In contrast. 



iNissen et al.l (2003a) provide strong evidence for a con- 
stant [S/Fe]~ 0.35 for stars with [Fc/H] < — 1 , and a 
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quasi-linear rise to this value from [Fe/H] = 0. This dis- 
crepancy among the observational results for [Fe/H] < — 1 
deserves to be resolved before S abundances are added 
to the observational inventory relevant to interpretations 
of Galactic chemical evolution. This paper describes our 
attempt to resolve the discrepancy. 

A problem with obtaining the sulfur abundance of 
metal-poor stars is the paucity of suitable atomic lines. 
Candidate lines of Si lie in the near- infrared (NIR): 
multiplets of high-excitation lines occur near 8694 A 
and at 9212 - 9238 A. The drawback with the Hues at 
8694.0 A and 8694 . 6 A is the ir weakness in halo stars . 
llsraelian fc Rebolol l|200j) and lTakada-Hidai et all l|2002ll 
observed these lines. On the other hand, a disadvantage 
with the much stronger triplet at 9212.9, 9228.1, and 
9237. 5 A is the he avy interference by the telluric H2O 
lines. iNissen et al.l (,2003a) observed these AA9213 — 38 
lines and corrected for the H2O line absorption. 

In this project, we chose to observe the AA9213 — 38 
lin es in a number of th e stars observed at AA8694 — 95 A 
by iTakada-Hidai et al.l ll2002l i n the following shortened 
as TH) and llsraelian fc Rebolol l)200li abbreviate as IfcR) . 
Our stars sample the metallicity range down to [Fe/H] ~ 
—3, and consist of thr ee from the sam ple of IfcR, three 
from TH, and four from lFrancoisI l)l988j) . whose equivalent 
widths were reanalyzed by TH and IfcR. Our primary goal 
was to check the high [S /Fe] reported previously from the 
AA8694— 95 lines. This check is executed using high quality 
spectra of the AA9213 — 38 hues analyzed using model 
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Fig. 1. In the upper panel, the full line displays the 
continuum-normalized spectrum of our giant star HD 2665 
(including telluric lines, mostly water vapor lines) and, as 
dots, the continuum-normalized spectrum of our rapidly 
rotating B-star, which shows only telluric absorption lines. 
In the lower panel, the HD 2665 spectrum is divided by 
this telluric spectrum, resulting in a pure stellar spectrum 
of HD 2665. The sulfur lines and the Paschen C hydrogen 
line emerge clearly after the division. Magnesium and iron 
lines are also marked. The abscissa and ordinate scales are 
the same in both panels. Observe that the local continuum 
of the HD 2665 spectrum is fitted with a low order func- 
tion, to avoid eliminating the wings of the hydrogen line. 
This will be even more important for the dwarfs 



atmospheres computed for the same effective temperature 
and surface gravity as used in the pioneering studies from 
which our star-list is made. As a secondary goal, we sought 
to gain an indication of how [S/Fe] behaves among halo 
stars of different metallicities. 

2. Observations 

The investigated stars, all lying in the northern sky, are 
presented in Table ^ I&R determined the sulfur abun- 
dance for 8 new stars, including two upper limits, with 
the 4.2 m William Herscliel Telescope on La Palma. We 
have analyzed three of these, namely HD 2665, HD 19445, 
and HD 201889. TH made new determinations of the sul- 
fur abundance in 6 stars (with one upper limit) using the 
HIRES spectrograph on the Keck telescope on Hawaii. Of 
these we have observed three (HD 84937, HD 88609, and 
HD 165195). Furthermore, we have observed HD 111721, 
HD 94028, HD 132475, and HD 201891, which were origi- 
nally observed and analyzed bv lFrancoisI l)l988|) . and then 
reanalyzed both by TH and I&R. 

The stars, which span the metallicity range of —2.75 
< [Fe/H]-'^< —0.7, were observed in observing runs from 
November 29 and 30, 2001 and April 27, 2002 at the 
W. J. McDonald Observatory using the 2.7-m Harlan J. 
Smith telescope with the high -resolution 2dcoude cross- 
dispersed echelle spectrograph ijTull et all [19951. The re- 
solving power was approximately R — 60 000, a value de- 
termined from the fwhm of thorium lines in the wave- 
length calibration frames. Full spectral coverage from 3400 
to 10 900 A can be obtained in two exposures. However, 
we observed every star in only one setting of the echelle, 
thus retrieving the entire wavelength range but with gaps. 
Our bandpass includes the sulfur AA9213 — 38 lines but 
not the AA8694 - 95 fines. We used the TK3 detector, a 
2048 X 2048 pixel^ CCD with 24/zm pixels. The observing 
times, ranging from a half-an-hour to 3 hours for the F, 
G, and K stars, are given in Table To avoid a large 

^ The bracket notation is defined as follows: 
[A/B] = log(A/B), - log(A/B)0 



number of cosmic hits, individual exposures were limited 
to a maximum of 30 minutes. 

The observed CCD data was processed with the reduc- 
tion package IRAF to retrieve one-dimensional, continuum 
normalized, and wavelength calibrated, pure stellar spec- 
tra. In the wavelength calibration a root-mean-square of 
the fit of less than 4 mA is achieved. Multiple exposures 
of the same star were added in the form of 1-dimensional 
spectra by weighting the different spectra by their mean 
count levels. When several frames are available a cosmic- 
ray-hit rejection algorithm was applied. 

The local continua of the F, G, and K star spectra were 
fitted and normalized by a low-order Legendre function in 
order not to remove the wings of the stellar Paschen line. 
As can be seen in Figure one of the sulfur lines lies in 
the wing of this hydrogen line. The hydrogen lines in the 
dwarf spectra have wings stretching over a considerable 
wavelength range and we therefore have to proceed with 
caution when fitting the continuum. 

The 9200 A region of interest here is marred by nu- 
merous telluric absorption lines (Figure^. To eliminate 
the telluric contribution to the spectra, we observed fast- 
rotating-B-stars with vsim = 150 — 350 kms^^. The in- 
trinsic, B-star spectra are featureless except for a more 
or less pronounced, rotationally broadened Paschen C H i 
line at 9229.02 A, which means that they can be used as 
a template for the telluric absorption lines. For these B- 
star spectra, we used a spline function of high order to fit 
and normalize the spectra, but also to eliminate the broad 
Paschen C Hi feature and occurring fringes. The hot stars 
were observed close in time and at a similar airmass as 
the programme stars, and were required to have a signal- 
to-noise ratio of several hundred. To achieve this in the 
9200 A region, exposure times from tens of seconds up to 
an hour were required, depending on the brightness of the 
individual star and the specific spectral type. The reduced 
B-star spectra achieved in this manner are assumed to con- 
tain purely telluric absorption lines. The programme spec- 
tra were subsequently divided by a scaled version of the 
'telluric' spectra. The divided, reduced spectra are shown 
in Figures El and 13 As can be seen in these Figures, the 
telluric lines are astonishingly well eliminated while the 
stellar sulfur lines show up clearly. While our procedure 
of eliminating the telluric lines adds to the noise, the re- 
sultant signal-to-noise ratio of the ratioed stellar-spectra 
is satisfactory for our purpose. 

The equivalent widths of the sulfur lines at 9212.9 A 
and 9237.5 A, which are considered to be unblended, were 
measured by Gaussian fits or by total integration of the 
line. The equivalent width of the 9228.1 A line is not mea- 
sured (except for HD 111721) since it lies in the stellar 
Paschen C Hi line. The measured equivalent widths of the 
two (three) sulfur lines, the equivalent widths of the Si 
8694.6 A line supplied in the literature (I&R; TH), and 
our measured equivalent widths of six Fe ii lines (lying at 
4416.8, 5264.8, 5284.1, 6149.3, 6432.7, and 6456.4A) are 
presented in Table [21 We note that the measured equiv- 
alent widths of the sulfur lines at 9212.9, (9228.1), and 



N. Ryde and D. L. Lambert: On the galactic chemical evolution of sulfur 
Table 1. Investigated stars 



3 



Star Spectral V Tea log((?) [F^/H] ^micro ^macro Date of Exposure time 

type^ [^-oiff] [K] (cgs) [km s^^] [km s^"""] observations [s] 

Giants and subgiants 

HD 2665 G5III 7.8 4990 ± 70 2.50 ±0.2 -1.74 ± 0.15 1.5 ± 0.5 5.5 ± 0.5 Nov. 29, 2001 7,200 

HD 88609 G5IIIw 8.6 4570 ± 100 0.75 ± 0.15 -2.85 ± 0.15 1.9 ± 0.5 7.0 ± 0.5 Apr. 27, 2002 3,600 

HD 111721 G6V 8.0 5010 ± 100 2.31 ± 0.15 -1.27 ± 0.15 1.2 ± 0.5 4.5 ± 0.5 Apr. 27, 2002 3,600 

HD 165195 K3p 7.3 4190 ± 100 1.00 ±0.15 -1.75 ±0.15 1.3 ± 0.5 7.5 ± 0.5 Apr. 27, 2002 1,800 

Dwarfs 

HD 19445 sdG5 8.0 5810 ± 150 4.46 ± 0.2 -1.90 ± 0.15 1.5 ± 0.5 4.5 ± 0.5 Nov. 29, 2001 10,800 

HD 84937 sdF5 8.3 6300 ± 100 3.97 ±0.15 -2.06 ±0.15 1.1 ±0.5 5.0 ±0.5 Nov. 30, 2001 9,000 

HD 94028 F4V 8.2 5980 ± 100 4.30 ± 0.15 -1.35 ± 0.15 1.5 ± 0.5 4.5 ± 0.5 Apr. 27, 2002 3,600 

HD 132475 F5/F6V 8.6 5810 ± 100 3.91 ± 0.15 -1.38 ± 0.15 1.9 ± 0.5 4.5 ± 0.5 Apr. 27, 2002 3,600 

HD 201891 F8V- VI 7.4 5880 ± 100 4.25 ± 0.15 -1.03 ± 0.15 1.5 ± 0.5 4.0 ± 0.5 Apr. 27, 2002 1,800 

HD 201889 GIV 8.0 5615 ± 100 4.24 ±0.2 -0.71 ± 0.15 1.5 ± 0.5 4.0 ± 0.5 Apr. 27, 2002 1,800 

From Simbad (http://simbad.u-strasbg.fr) 



9237.5 A range from approximately 10 to 100 mA, which 
is suitable for a reUable abundance analysis. Observe also, 
that the equivalent widths are an order of magnitude 
larger than those determined for the line at 8694.6 A. The 
8694.0 A line is a factor of five weaker still and unmeasur- 
able in all the metal-poor stars we have been investigating, 
except HD 201889. 

We note that the TK3 detector used consists of a thick 
chip, which is not severely affected by fringing. The exist- 
ing fringing cancels well using a flat field - the top panel is 
Figure n compares two such flat- fielded spectra. Fringing 
is eliminated to a high degree by using a hot star observed 
at the same time and the lower panel of Figure ^ shows 
just such a spectrum. 

In Figures 121 and 131 there exist low points which orig- 
inate from the cancellation of the very strongest tel- 
luric lines. None of these affect the SI lines. For exam- 
ple, in the spectrum of HD 19445 in Figure El only the 
9212.9 A line lies in a wing of a modestly strong tel- 
luric line. The other ones are free. The abundance de- 
rived from the free 9237.5 A line is log As = 5.51 (the 
logarithmic number-abundance of sulfur relative to hy- 
drogen, log As = logiVs/iVn, normahzed on a scale where 
log Ah = 12.00), whereas the ratioed 9212.9 Aline gives 
5.49. The synthetic spectrum which fits these two lines, as 
well as the 9228.1 A line is calculated with an abundance 
of 5.50. The fact that the three sulfur lines are synthesised 
well with one sulfur abundance shows that the elimina- 
tion of the telluric lines produces a reliable result, clearly 
within the other uncertainties. We estimate the error in 
the equivalent widths to be less than 5%, implying an un- 
certainty in the derived sulfur abundance of the order of 
0.05 dex, as estimated from the individual measurements. 

In some regions, for example around 9207 A in the 
spectrum of HD 111721 in Figure El and around 9245 A 
in the spectrum of HD 132475 in Figure the elimina- 



tion of the telluric lines is hampered due to closely lying 
telluric lines with wings blending with each other. The 
equivalent widths of the sulfur lines are, however, mea- 
sured by drawing a local continuum which should reduce 
systematic errors due to this effect. 

3. Analysis 

We analyze our data by modelling the stellar atmosphere 
and requiring that the measured equivalent widths are re- 
produced for the two unblended S i A9213 and A9238 lines. 
The mean sulfur abundance yielded in this way is sub- 
sequently used in calculating a synthetic spectrum for a 
given atmosphere of the entire region (9200-9500 A). The 
synthetic spectrum is thereafter convolved with a macro- 
turbulence function in order to fit the shapes and widths of 
the lines, see Figures|21and|31 In this section we will discuss 
the model atmospheres, the stellar parameters including 
uncertainties, the line data, and the spectrum synthesis. 

3.1. Model atmospheres 

I&R and TH used ATLAS9 model atmospheres whereas 
we use model atmospheres provided by the MARCS code. 
The M ARCS code was first developed bv lGustafsson et alJ 
and has been successively updated ever since. 
These hydrostatic, plane-parallel model photospheres are 
computed on the assumptions of Local Thermodynamic 
Equilibrium (LTE), homogeneity and the conservation of 
the total flux (radiative plus convcctivc; the convective 
flux being computed using the mixing length formulation) . 
Data on absorption by atomic species are c ollected from 
the VALD database d Piskunov et alJ . Il995f) and Kurucz 
(1995, private communication). Absorption by molecules 
is included but quite unimportant for our stars. 

We have computed model atmospheres in plane- 
parallel geometry which is an excellent approximation for 
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Table 2. Measured equivalent widths given in inA 



Wsi Ws, Ws: Wpeii Wpen Wpen Wpen Wpen ^Feii 

9212.9 A 9228.1 A 9237.5 A 8694.6 A 4416.8 A 5264.8 A 5284.1 A 6149.3 A 6432.7 A 6456.4 A 
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From Israelian & Rebolo (2001) 
^ From Takada-Hidai et al. (2002) 
From Frangois (1988) 



Fig. 2. Continuum-normalized spectra of our programme giants: HD 2665, HD 88609, HD 111721, and HD 165195. 
Telluric absorption lines are ratioed out, but some residual signatures are visible. In the spectra of HD 2665 and 
HD 111721 several Fei lines and a Mgii line are discernable and in the spectrum of HD 165195 the most conspicuous 
Pel line is detected. In the lower right panel in Figure |31 identifications of all observable metal lines in this wavelength 
region can be found. Our best model spectra are also plotted with full lines. See text for a comment on the bad fit 
of the hydrogen Paschen line in the spectrum of HD 165195. After the names the star-model's temperature, \ogg, 
metallicity, and microturbulence are indicated 

Fig. 3. Continuum normalized spectra of our programme dwarfs: HD 19445, HD 84937, HD 94028, HD 132475, 
HD 201891, and HD 201889. The telluric absorption lines are removed as well as possible, but some residuals remain, 
most prominently in the spectrum of HD 19445 and HD 84937. Metal lines have been identified in the spectra of the 
four lowest panels. In the spectrum of the most metal-rich star, HD 201889, all the identified Pe i, Mg ii, and S I 
lines are marked. The full lines represent our modelled spectra. After the names the star-model's temperature, logg, 
metallicity, and microturbulence are indicated 



the dwarfs. To test the validity of the approximation for 
the giant stars we calculated a model atmosphere in spher- 
ical geometry for the giant HD 88609 and compared it 
with a plane-parallel model. The equivalent widths of the 
sulfur lines are changed by only a few percent. We have, 
therefore, calculated model atmospheres for all our stars 
in plane-parallel geometry. 

In our model atmospheres the enhancement of C, O, 
Ne, Mg, Si, S, Ar, Ca, and Ti (that is the a elements) 
was assumed to be [a/Pe] = 0.4 except for HD 201889 for 
which it was assumed to be [a/Pe] = 0.2, due to its higher 
metallicity. The derived S abundance is quite insensitive 
to the a enhancement. 

Using their equivalent widths, line parameters and set 
of fundamental stellar parameters, we are able to repro- 
duce the derived abundances in I&R and TH. Our use 
of MARCS model atmospheres lowers the abundances by a 
few hundredths of a dex relative to those obtained from 



ATLAS9 models used previously. Clearly, the choice of AT- 
LAS9 vs. MARCS is not an important factor affecting the 
derived abundances. 

3.2. Stellar parameters 

We have chosen the same stellar, fundamental parame- 
ters for the model atmospheres (actually Tcff, ^og{g), and 
Cmicro) SiS thosc determined in I&R and TH. However, the 
metallicities of the stars ([Fe/H]) we redetermined from 
six singly ionized iron lines. The fundamental stellar pa- 
rameters including our determined metallicities and are 
given in Tabled 

As was noted above, the stars HD 111721, HD 94028, 
HD 132475, and HD 201891, which w ere originally ob- 
served and analyzed bv lPrancoisI l)l988() . have been reana- 
lyzed both by TH and I&R. For the determination of the 
effective temperatures and surface gravity, TH took into 
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account interstellar reddening which resulted in noticeable 
adjustments of the parai neters. We will u se parameters de- 
termined by TH for the iFrancoisI l)l988l) stars. 

TH determined the effective temperatures of 
HD 84937, HD 88609 HP 165195, and HD 111721 
from the lAlonso et all l|l999a|) calibration, based on the 
infrared flux method (irfm), also taking the interstellar 
reddening into account. For the reanalysis of the dwarfs 
HD 94028, HD 132475, and HD 201891, TH used the 
empirical temperatu re scale for dwarfs as formulated by 
lAlonso et al"l l|l996|) . Our new metallicities do not alter 
this temperature determination significantly. 

Furthermore, TH calculated the surface gravity us- 
ing Hipparcos parallaxes. The microturbulence was deter- 
mined for HD 84937, HD 88609, and HD 165195 from 
Fe I lines . Microturbulence for HD 111721 is taken from 
iR.van fc Lambert (1995) . and for the dwarfs HD 94028, 
HD 132475, and HD 201891 it was calculated f rom t he 
empirical formula provided by lEdvardsson et af] l)l993|) . 

For the giant star HD 2665, the effective temperature 
was det ermined by I&R based on the irfm in lAlonso et alJ 
lll999bt) and the surface gravity (log 5) from a non-LTE 
study of iron (for references see I&R). For HD 19445 
and HD 201889 the temperature is based on the irfm 
raelian et al. UI99?), and log q from an non-LTE analy- 
sis of iron bv lThev cnin fc Idiart (1999). The microturbu- 
lence we have assumed to be Xmicro = (1.5 ± 0.5) km s^^ 
for these three stars. 



3.3. Line data 

The high-excitation {x = 6.5 eV) S i lines used in this 
paper are due to the multiplet As^S^ — Ap^P, while 
the lines at 8694.0 and 8694.6 A are from the multiplet 
Ap^P — Ad^D^ . The line parameters (wavelength, excita- 
tion energy, line strength) of our S I lines and six Fe 11 lines 
are provided in Table |3| 

In the following, we h ave used a solar s ulfur abun- 
dance of log £0(8) = 7.20 l|Chen et all l2002|) as a refer- 
ence point. 



Table 3. Line parameters 



3.4. Spectral synthesis 

The sulfur abundance is originally found by reproducing 
the equivalent widths of the unblended sulfur A9213 and 
A9238 lines. This is done with the radiative transfer rou- 
tines of the MARCS codes and yields a mean sulfur abun- 
dance with a standard deviation of 0.01 to 0.1 dex. (This 
gives an assessment of the random measuring uncertain- 
ties.) Synthetic spectra around 9200 — 9250 A are pro- 
duced by computations of the radiative transfer through 
the model atmospheres, using the mean sulfur abundance 
and using our line data for all three sulfur lines and the 
Paschen C Hi line. We calculate the radiative transfer for 
points in the spectrum separated by AA ^ 0.6 kms^^ (cor- 
responding to a resolution of A/AA ^ 460 000). We sub- 
sequently convolve the spectra with a macro-broadening 



Wavelength 


Ecxc 


log gf 
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(eV) 


(cgs) 


S I LINES 








6.525 


0.43^ 




6.525 


0.25^ 




6.525 


0.04^ 


Fe 11 LINES 






4416.830 


2.778 


-2.80'' 


5264.812 


3.230 


-3.19'' 


5284.109 


2.891 


-3.29'' 


6149.258 


3.889 


-2.82'' 


6432.680 


2.891 


-3.67'= 


6456.383 


3.903 


-2.25" 



Line data taken from the nist database 
'' From Thoren, Edvardsson, & Gustafsson, in preparation 
" From Thoren & Feltzing (2000) 



function, assuming a Doppler-shift distribution for both 
radial and tangential velocity components as specified 
by the 'radial- tangential' m odel for the macroturbulence, 
for details see iGravl lll992t) . The FWHM velocities for the 
macroturbulence broadening are given in Tableland were 
derived by requiring that the shape and widths of the sul- 
fur lines should match the observations. 

Note, that the sulfur A9228 line, which lies in the 
Paschen C, Hi line, will also be calculated (including the hy- 
drogen opacity) by the synthetic spectrum program. For 
a few of the dwarfs the wings of the Paschen Q Hi line 
also interfere, to various extents, with the two sulfur lines 
which were assumed to be unblended. This means that the 
measured equivalent widths could be wrong in these cases, 
underestimating the sulfur abundance. Therefore, the sul- 
fur abundances were adjusted in order that the synthetic 
spectrum, including the hydrogen-line wings, should fit 
the observations. The size of these changes were, however, 
not more than approximately 0.05 dex in the sulfur abun- 
dance. 

In Figures 13 and |21 we have plotted the final syn- 
thetic spectra together with our observations. The mod- 
elled spectra are shifted according to the observed radial 
velocities of the stars. From the Figures, we see that all 
sulfur lines are synthesized convincingly. Thus, the infor- 
mation from all three sulfur lines has been used in deter- 
mining the sulfur abundance. The synthetic spectra only 
take the sulfur lines and the hydrogen line into account. 
The other metal lines in the region were not included in 
the calculation. 

The calculations of the wings of the hydrogen line yield 
a satisfactory match to the observations, which can be 
seen, in particular, in the dwarf spectra (see Figure 13). 
Recall that the continuum normalization was made with 
a low order function only, in order to fit the continuum 
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Table 4. Effects on logarithmic abundances derived when 
changing the fundamental parameters of the model atmo- 
spheres. Two stars of different parameters are presented. 
The parameters (Toff/ log g/[Fe/H]/^t) are, respectively: 
4190/1.0/ - 1.75/1.3) and 6300/3.97/ - 2.06/1.1) 



Table 5. Abundances of Fe ii and S i, the latter based on 
the AA9213 - 38 lines 





UncGrtctinty 


A fc; /wi 


ZA [o / r cj 
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— U.i 1 
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HD 84937 


ST^g = +100 K 


-0.03 


-0.04 




51og5 = +0.2(cgs) 


+0.06 


-0.01 




5[Fe/H] = +0.15 


±0.00 


±0.00 




^Cmicro = +0.5 


-0.06 


-0.05 



Star 


log^Fell 


[TT^TT /VI] 

[r eil / jlj 


log As 


P/flJ 


re 


Giants and subgiants 










HD 2665 


5.76 


— 1.74: 


5.60 


— 1.60 


0.14 


HD 88609 


4.65 


—2.85 


4.74 


—2.46 


0.39 


HD 111721 


6.23 


-1.27 


6.25 


-0.95 


0.32 


HD 165195 


5.75 


-1.75 


5.97 


-1.23 


0.52 


Dwarfs 












HD 19445 


5.60 


-1.90 


5.50 


-1.70 


0.20 


HD 84937 


5.44 


-2.06 


5.27 


-1.93 


0.13 


HD 94028 


6.15 


-1.35 


5.93 


-1.27 


0.08 


HD 132475 


6.12 


-1.38 


5.97 


-1.23 


0.15 


HD 201891 


6.47 


-1.03 


6.26 


-0.94 


0.09 


HD 201889 


6.79 


-0.71 


6.69 


-0.51 


0.20 



over a limited range, so spurious features (such as possible 
residual fringing) are not taken out. The hydrogen line 
cores are expected to be subject to departures from LTE, 
which could lead to poor fits. The unsatisfactory fit of the 
hydrogen-line core in the giant HD 165195, for example, 
could be improved but would require higher temperatures, 
of the order of 400 K. 

3.5. Uncertainties in the stellar parameters 

The uncertainties in the fundamental parameters are 
adopted from their sources: for HD 2665, HD 19445, and 
HD 201889 from I&R; for HD 84937, HD 88609, and 
HD 165195 from TH; and for HD 111721, HD 94028, 
HD 132475, an d HD 201891 from the reanalysis of the 
iFrancoisI (|l988|) data by TH, see Table □ For HD 2665, 
HD 19445 and HD 201889, we assess the uncertainties 
in the microturbulence to be of the same order as esti- 
mated in TH. Furthermore, we judge the uncertainties 
in the macroturbulence to be of the same order, that is 
±0.5kms~^. The uncertainties in our new determination 
of the metallicity is assessed to be less than ±0.15 dex. 

The propagation of these uncertainties into the deter- 
mination of the sulfur abundance is calculated by chang- 
ing the fundamental parameters of the model atmosphere 
with the estimated uncertainties and then running the 
synthetic spectrum program to obtain the change in sul- 
fur abundance for a given equivalent width. Table 01 shows 
the consequence of the uncertainties in the effective tem- 
perature, surface gravity, metallicity, and microturbulence 
on the [S/H] and [S/Fe] ratios for one giant and one 
dwarf. For the coolest stars (e.g. HD 165195), the domi- 
nant source of uncertainty is the 100 K estimate for the 
temperature error. The uncertainties in the sulfur abun- 
dance, due to temperature uncertainties, are of the order 
of 0.04 dex or less for the dwarfs, and for the giants, typ- 
ically 0.1 dex. 



The total uncertainty in the [S/Fe] due to uncertainties 
in the fundamental parameters is, in general, of the order 
of 0.05, but could be as high as 0.08 in specia l cases. This 
is the same conclusion as lNissen et al.l l|2003a^ arrive at in 
their investigation of the sulfur abundance in solar-type 
stars. 

3.6. Sulfur and iron abundances 

The resulting sulfur abundances for our stars are presented 
in Table|5las log Ag, the logarithmic abundance relative to 
the solar value of log As,© — 7.20, i.e. [S/H], and the log- 
arithmic abundance relative to the solar value normalized 
to the metallicity ([S/Fe]). Also, displayed in the Table 
are the iron abundances calculated from the measured 
equivalent widths of Fe ii lines shown in Table [21 The to- 
tal uncertainty, including effects of statistical measuring- 
uncertainties and uncertainties in the model parameters, 
is estimated to be on the order of ±0.15 dex. 

In their non-LTE analysis of neutral sulfur in envi- 
ronments relevant also in our study, TH find inconse- 
quential effects for the AA8694 — 95 lines, ranging from 
0.00 to 0.08 dex in the sulfur abundance. iNissen et alJ 
l(2003tJ ^ find in their study of the sulfur abundance for 
similar stars, that the sulfur abundance they derive from 
the AA8694 — 95 lines match the abundance derived from 
the AA9213 — 38 lines very well, indicating that depar- 
tures from LTE should be small also for the latter lines. 
It should be noted that neutral sulfur atoms represent the 
main ionization state throughout most of the atmosphere, 
and most importantly in the line-forming regions. Thus, 
we judge that the NLTE effects should be small, iusti- 
fying a standard LTE analysis. Furthermore. iNissen et all 
(2003^) investigate and find that the 3D effects^ for dwarfs 
will not severely alter the [S/Fe] ratio (less than 0.05 dex). 

^ By '3D effects' is meant differences in, for example, derived 
abundances when using atmospheric models including 3D, hy- 
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Table 6. Sulfur abundances compared to those from the 
literature which are based on the AA8694 — 95 lines 



Star 


log As 


log j4s 






AA9213 - 38 lines 


AA8fi94 — 95 lines 






this work 


literature 




Giants and subgiants 






HD 2665 


5.60 


5.89 


(1) 


HD 88609 


4.74 


< 5.49 


(2) 


HD 111721 


6.25 


6.3/6.02='/-'' 


(3) 


HD 165195 


5.97 


5.98 


(2) 


Dwarfs 








HD 19445 


5.50 


5.97 


(1) 


HD 84937 


5.27 


5.70 


(2) 


HD 94028 


5.93 


6.3/6.06^/6.09'' 


(3) 


HD 132475 


5.97 


6.4/6.00^/6.08'' 


(3) 


HD 201891 


6.26 


6.6/6.29^/6.33'' 


(3) 


HD 201889 


6.69 


6.83 


(1) 



ril: Israelian fc Rebolo ('2001') 

(2) : Takada-Hid ai et al. (2002) 

(3) : Frangois (1988) 

^ A reanalysis by Takada-Hidai et al. (2002) using Frangois's 

(1988) equivalent widths 
'' The corresponding reanalysis by Israelian & Rebolo (2001) 



Fig. 4. This figure shows the Galactic chemical evolu- 
tion of sulfur. The iron and sulfur abundances, as we 
have determined it from analyzing Fe ii lines and sulfur 
lines in the NIR, are presented with star symbols. Typical 
errorbars for our measurements of ±0.15 dex, both for 
[S/Fe] and [Fe/H], are indicated. Circles indicate the sulfur 
and iron abundances derived byti rraeli an_ fc Rebo lo ( 2001, 
and sq uares display values found by iTakada-Hidai et 
l)2002() (including an upper limit which is illustrated with 
a triangle). For the 6 s tar we have in common wit h 
llsraeUan fc Rebolol l|200lh and lTakada-Hidai et alJl(20n2ll . 
we have connected our new determinations with their mea- 
surements, which are based on the sulfur AA8694 — 95 
lines. We have also p lott ed sulfur abundanc es derived by 
iNissen et aP ()200 3a') and lChen et all l)2002l) for southern 
stars. These are plotted with crosses 



Therefore, we will neither make any corrections for non- 
LTE effects nor for 3D effects. Our [S/Fe] values are plot- 
ted versus [Feii/H] in Figure^ It is clear that departures 
from LTE and the 3D effects have very similar (small) ef- 
fects on the abundance derived from the AA9213 - 38 lines 
and that from the AA8694 — 95 lines. 



drodynamic modelling of the convection compared with tradi- 
tional atmospheric models, which describe convection by the 
mixing-length approximation. 



Since sulfur abundances have been obtained previously 
for each of our stars using the AA8694 — 95 lines, it is of 
considerable interest to compare S (and Fe) abundances. 
This comparison is made in Tables IHl and 

Three stars - one giant and two dwarfs - are in com- 
mon with I&R. Our sulfur abundances are systematically 
smaller: A log As — log A5 (ours) — log Ag (theirs) ~ —0.29 
(HD 2665), -0.47 (HD 19445), and -0.14 (HD 201889). 
Recall that we have used the same stellar parameters. It 
is likely to be significant that the agreement is best for 
the star (HD 201889) for which the AA8694 - 95 Hues are 
strongest. The difference is in fact —0.14 dex also when the 
abundance of I&R is corrected to the nist gf- value (abun- 
dance increased by 0.02 dex) and to the MARCS model 
(abundance decreased by 0.02 dex). We suppose that the 
larger differences for the two more metal-poor stars re- 
fiect the greater uncertainty in measuring very weak S i 
A8695 line (equivalent widths of about 3 mA) . Alternative 
choices of model atmosphere parameters cannot erase the 
abundance difference resulting from the choice of S i lines. 
The corresponding As for [Fe/H] are +0.26 (HD 2665), 
-0.02 (HD 19445), and +0.23 (HD 201889). I&R indicate 
that the iron abundance determination is based on Fei 
lines and a non-LTE analysis. We greatly prefer to derive 
the Fe abundance from Fell lines, for which non-LTE ef- 
fects are slight, than to use Fei lines and an uncertain 
non-LTE correction. Fe 11 lines are, however, more sensi- 
tive than Fe i lines to uncertainties in the surface gravity, 
but this is also true for S i lines compared to S 11. Thus, 
this sensitivity should cancel out in the Si/Feii ratio. 

There are large differences in the [S /Fe] ratios obtained 
by us and I&R: 0.14 vs. 0.69 (HD 2665), 0.20 vs. 0.65 
(HD 19445), and 0.20 vs. 0.57 (HD 201889) arising from a 
higher S abundance in each case and a lower Fe abundance 
in two cases. In short, we do not confirm the high [S/Fe] 
proposed by I&R. 

Three of our stars were observed and analyzed pre- 
viously by TH. For HD 88609, the most metal-poor star 
([Fe/H] = —2.8) of the trio, their upper limit to the S 
abundance is a safe 0.7 dex higher than our measured 
abundance. For HD 84937 with [Fe/H] = -2.1, the sul- 
fur abundance derived by TH is 0.45 dex larger than ours 
but this difference is probably due to the uncertain equiv- 
alent width of the A8695 line; TH note that their spec- 
trum is of low signal-to-noisc ratio and marred by resid- 
ual fringing. In the case of HD 165195, our and their S 
abundances agree well. These comparisons use LTE abun- 
dances derived by TH; the corrections for non-LTE effects 
are small. The LTE Fe abundances derived by TH from a 
set of Fell lines are in good agreement with ours: As are 
-0.12, -1-0.12, and +0.05. The [S/Fe] values for the two 
stars for which TH detected the A8695 line in the sense 
ours vs. theirs, are 0.52 vs. 0.64 (HD 165195) and 0.13 vs. 
0.60 (HD 84937). TH detected the AA8694 - 95 fines in 
three stars not observed by us. Their detections appear 
quite secure. We note that the [S/Fe] ratio inferred by TH 
for these stars, with [Fe/H] ~ —1.4, average about 0.2 and 
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Table 7. The iron and sulfur abundance ratios ([Fe/H], [S/H], and [S/Fe]) from our work and from the literature for 
the AA8694 - 95 hues 



Star [Fe/H] [S/H] [S/Fe] [Fe/H] [S/H] [S/Fe] reference 



based on the AA9213 — 38 lines study based on the AA8694 — 95 lines 



Giants and subgiants 



HD 2665 


-1.74 


-1.60 


0.14 


-2.0 


-1.31 ±0.14 


0.69 


(1) 


HD 88609 


-2.85 


-2.46 


0.39 


-2.73 


< -1.72 


< 0.99 


(2) 


HD 111721 


-1.27 


-0.95 


0.32 


-1.57/-1.29''/- 


-0.86/-1.19^/- 


0.71/0.07 ± 0.41='/- 


(3) 


HD 165195 


-1.75 


-1.23 


0.52 


-1.87 


-1.23 


0.64 ±0.09 


(2) 


Dwarfs 
















HD 19445 


-1.90 


-1.70 


0.20 


-1.88 


-1.23 ±0.12 


0.65 


(1) 


HD 84937 


-2.06 


-1.93 


0.13 


-2.11 


-1.51 


0.60 ±0.06 


(2) 


HD 94028 


-1.35 


-1.27 


0.08 


-1.48/-1.26''/ - 1.31*^ 


-0.91/-1.15''/ - 1.12*^ 


0.57/0. ll'^/O.ig'' 


(3) 


HD 132475 


-1.38 


-1.23 


0.15 


-1.37/-1.23''/ - 1.37*^ 


-0.81/-1.21''/ - 1.13*^ 


0.56/O.Ol='/0.24'' 


(3) 


HD 201891 


-1.03 


-0.94 


0.09 


-1.25/-0.98''/ - 0.87*^ 


-0.61/-O.92''/ - 0.88*^ 


0.64/0. 05="/ - O-Ol*" 


(3) 


HD 201889 


-0.71 


-0.51 


0.20 


-0.94 


-0.37 ±0.1 


0.57 


(1) 



fl): Israelian fc Rebolo (20011 

(2) : Takada-Hidai et al. (2002') 

(3) : Frangois (1988) 



A reanalysis by Takada-Hidai et al. (2002) using Frangois's (1988) equivalent widths 
^ The corresponding reanalysis by Israelian & Rebolo (2001) 



fit well with our suggested run of [S/Fe] with [Fe/H], see 

Figure 01 

We chose four stars from ' Francois! ijlQSSl) . This quartet 
was reanalysed by TH and IfcR using equi valent widths of 
the A8695 line reported bv iFrancoisI l) 1 QSSj) . TablelHlshows 
that, for the dwarfs, these reanalyses and our abundance 
from the NIR lines are in good agreement. The agreement 
extends to [S/Fe], cf. Table [7| The typical total uncertain- 
ties quoted by TH, for their dwarfs, for the [S/Fe] ratio 
is 0.15 dex. However, for their analysis of HD 111721 the 
uncertainty is larger (0.4 dex). 



4. Discussion 

Our observations show that the Si NIR lines are de- 
tectable and S abundances determinable to [Fe/H ] —3 , 
and, thus, confirm the conclusions of Nis sen et al.l l)2003atl 
from their vlt/uves observations of the same lines in 
southern dwarfs and subgiant s. In Figure IH we show our 
results with those of .Nissen et al.l ll2003a^ for halo stars, 
also based on t he AA9213- 3 8 line s, and lChen etlH l|2002() 
for disk stars. IChen et all (|20P2|) used the AA8694 - 95 
lines and also lines at 6046, 6053, and 6757 A, but for 
their chosen disk stars these lines provide readily mea- 
surable lines in high-resolution spectra. In Figure 0] our 
values for the stars in common with I&R and TH are con- 
nected with theirs to guide the eye. For the stars we have 
in common with TH, we have, for consistency, chosen their 
[Fcii/HJlte and [S/Feii]LTE values. The figure illustrates 



the comparisons between our and the AA8694 - 95-based 
abundances discussed in the previous section. 

There is no evidence in our data for the rise of [S/Fe] 
with decreasing [Fe/H] earlier proposed by I&R and TH. 
I&R claimed a linear trend such that [S/Fe] ~ 0.9 at 
[Fe/H] = -2, and TH put the slope lower so that [S/Fe] 
— 0.5 at [Fe/H] = —2. We have suggested that the lower 
[S/Fe] now being obtained are due to a combination of 
overestimates of the strength of the very weak AA8694 — 95 
lines in the most metal-poor stars, and of underestimates 
of the Fe abundances by I&R. We reiterate that the 
adopted effective temperatures and surface gravities are 
those used by I&R and TH. Our primary goal was to check 
the latter authors' S abundances, as directly as possible, 
by using AA9213 - 38 hues rather than the AA8694 - 95 
lines used by them. Definition of the run of [S/Fe] with 
[Fe/H] was a secondary goal whose full achievement de- 
mands additional stars. As can be seen in Figure 01 the 
evidence that [S/Fe] is nearly independent of [Fe/H] for 
[Fe/H] <— 1 is greatly s trengthened o n addi ng the sulfur 
abundances derived bv iNissen et alJ I 2fl03a|l . who find a 
constant value of [S/Fe] « 0.35 dex. 

There is a hint in Figure 01 that our [S /Fe] estimates 
may be sy stematically about 0.1 dex lower than those of 
iNissen et al.„ (2003a ) . This difference could be the result of 
an accumulation of minor differences in effective tempera- 
ture scales, surface gravity determinations, Fell gf-values, 
and others. We also note, that .Nissen et al.. (,2003al derive 
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a solar iron abundance 0.03 dex larger than the one we ap- 
ply. Considering this fact would diminish the discrepancy. 

Many studies of abundance of a elements in halo stars 
have shown that [a/Fe] at a fixed [Fe/H] is largely with- 
out co smic scatter. The S abundances of iNissen et all 
2003a^ confirm that this result may now be extended to 
sulfur. In our small sample, HD 165195 appears to de- 
part from the constant [S/Fe] indicated by other stars. 
This departure may possibly be due to adoption of incor- 
rect atmospheric parameters. We note that TH remark 
that the determination of the effective temperature for 
HD 165195 is uncertain. They consider a range of temper- 
atures (4131 — 4507 K) in which their adopted tempera- 
ture (and therefore the one used by us) is at the lower end 
of this range. Observe, that the sulfur abundances deter- 
mined for this star are the most sensitive to the effective 
temperature of our entire sample, which can be seen in 
Table 01 A temperature increase of 300 K would yield a 
[S/Fe] ratio close to our mean value. Therefore, our value 
might be seen as a determination which could be on the 
high side. We note, in passing, that an increase in tem- 
perature would also make the modelled hydrogen-line fit 
the observed line better. 

The exact run of the abundances with metallicity 
has consequences for theories and our understanding of 
Galactic chemical evolution and for the sites of forma- 
tion of the elements. So, what are the sites where sulfur 
was synthesized during the first few billion years of the 
Universe? a-capture elements are thought to be formed 
during explosions of supernovae (SNe) Type II. Owing to 
the short life-times of high-mass stars, the progenitors of 
SNe Type II, the abundances of elements synthesized by 
them and ejected into the interstellar medium, quickly 
reach a steady state and show a constant over-abundance 
in the ordinary [S/Fe] vs. [Fc/H] diagram. However, based 
on their data, TH and I&R discern a linear rise of the a- 
capture element sulfur for the halo phase of the Milky 
Way, thus showing a divergent behaviour compared to 
other a-capture elements, except possibly oxygen. If there 
is a rise in the abundances with decreasing metallicity, 
other processes for their formation have to be invoked. 
TH and I&R suggest a sulfur contribution from hyper- 
novae in the early galaxy as an explanation of the high 
values. Hypernovae are thought to be the explosion of an 
extremely massive star (several 100 M0). Another mecha- 
nism that could lead to a divergent evolutionary behaviour 
of oxygen and sulfur as compared to the other a-capture 
elements, is the fact that these two elements are volatile 
elements as opposed to the others which are refractory 
elements. This is important in the modelling of the trans- 
port and mixing of supernovae ejecta into the interstel- 
lar medium, with volatile elements experiencing a mix- 
ing time-scale, more than an order of magnitude faster. 
However, considering the results that .Nissen et al., (,,2003a} 
and we find, there does not seem to be a need for a dis- 
tinct evolutionary scenario for sulfur as compared to the 
other a. elements. 



5. Conclusions 

We have shown that reliable abundances of sulfur down to 
metallicities of [Fe/H] ~ —3 are obtainable using the S i 
lines lying at 9212.9, 9228.1, and 9237.5 A. In an attempt 
to resolve the discrepancy found in the literature, about 
the Galactic chemical evolution of sulfur for metallicities 
lower than [Fe/H]< —1, we have observed an ensemble 
of ten stars, both dwarfs and giants, previously analyzed 
with a different and weaker diagnostic; the S i lines at 
8694.0 and 8694.6 A. These previous investigations of the 
10 stars claimed a rise of the [S/Fe] ratio with decreas- 
ing metallicity. We are not able to confirm this rise. Our 
conclusion is instead that we corroborate the finding by 
iNissen et al.l l(2003al) indicating a similar behaviour of sul- 
fur to the other a elements. 

We suggest that the reasons for the difference in de- 
termined sulfur abundances and [S/Fe] are the strength 
of the lines analyzed and the way the metallicity is de- 
termined. First, we have used lines which are stronger by 
a factor of ten as compared to those used earlier. This is 
critical for metal-poor stars. Thus, we conclude that the 
smallest equivalent widths for the weaker A8695 line must 
have been overestimated in the previous studies, leading 
to an overestimation of the [S/H] ratio by up to 0.5 dex. 
Second, we have determined the metallicity from Fe 11 lines 
which should give a better determination of the metallic- 
ities than an non-LTE analysis of Fei lines. 

We also show that in spite of the numerous telluric, 
water-vapor lines in the 9200 A spectral region, a careful 
reduction of the data can provide clean spectra. We con- 
clude that the 9212.9, 9228.1, and 9237.5 A lines are the 
preferred ones to be used for abundance analyses of sulfur 
of halo stars. 
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